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CIE 1976 u-v f diagram 



u' 




CORM 2011 


NIST 2 


Duv often missing 


Lighting Facts Label 


Light Output 1 Lumens — 

tr^ hiflhtf tNi 
rumbef. Iba rrore Igtn is etnlLed 

Sttd'T'i « *+Kp**i Fiji ar 

Watts — 

AAmsufes BTKinJy iEquir*rJ to «t^i| 

tie pfcdutl The towesT Eta the 

istjwgy L-s*d 

Heoofad u Fnwr n^itaf ^ LM T- Tier. 

Lumens per — 

WatfcEffScacy 

MeesuRiB fifhceiH:^ The higbe^ Utb 

fumbv. if* nm«mc4Atihd pf&fcd 

aifrinw^Tit^ 1W : ura Lmi Tfo*-! 


I ESN A LH.79-200S 

WkfUTj MiMM HNP P itw* Nl 
^U0i^ll$ftifT^TTUi , V£ HUiHt d liDWJrtl 
jrd «!fH h rim Ibr ■ hf- mrpmtr 
d njgmSiw kujif 


RogialratiQii Number 
Model Number 
Type 


lighting facts 

" " A r*l3 UF iflTli Cl KH 


Liphi Omput (LumottsS 
Walls 

Luhlfrnft pet Willi; (EHiCdty) 


C«HPP ACCUPAGE 


PJr^mT K 


cu 



ST 


1 L>fl M r "‘a- 
1 i“j- »hiu. 

2fl00 (Warm Whrtfl) 


r 

| m-p^wra* 


cm*B“ 




itnx 


ft* 5 M 


iA* ^HOfl -.v 
i«« TlaU* 




VliM lupShnu-^Lni^ i v . Dfe. Li^- 7 ' flrfnmi GhrtJ* 


ijia -ll' l MCH1S«2tthP£>rri r THI' I 


Brand 


Color Rendering 
Index (CRI) 

r/vMSijnjs CAlCK MXlnSCy 

t&fr ImHir -1 T"> #rwi Id ^ L#n^i J i=$t. 

IpdOCU*! Drt-Sli COM IKWHU d CCfta 

Correlated Color 

Temperature (CCT) 

F.VjfliurflS Idjlil tttv 

"CMf M n#rt hfNi Mw WfWH (J6W-SH0 KJ 
#art f uj^jii r hi.j LflfT(iftiaiiiirt 1 7 r LC- 2C5 l£< K| 

C'Kf i>/v tap 4HK It >A\< Erin 

ztaVxri intrt" ir Mj4. . bd na r be i,'SF*3prak 

MHr ■gfaiiMi 


CCT and CRI do not tell the whole story of color quality 


CORM 2011 


^ ■ 


mm^m 


NIST 3 







CCT and CRI do not tell the whole story 




CIE 1 976 u’-v’ diagram 



u’ 




Triphosphor 
FL simulation 



Duv is another 
important 
dimension of 
chromaticity. 
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Duv defined in ANSI standard 


Closest distance from the Planckian locus on the ( u\ 2/3 v) 
diagram, with + sign for above and - sign for below the 
Planckian locus. (ANSI C78.377-2008) 


CIE 1960 u-v Chromaticity Diagram* 



Symbol: D uy 

CCT and Duv can specify 
the chromaticity of light 
sources just like (x, y). 

The two numbers (CCT, Duv) 
provides color information 
intuitively, (x, y) does not. 

Duv needs to be defined 
by CIE. 
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ANSI C78. 377-2008 Specifications for the 
chromaticity of SSL products 


CIE 1931 (x, y) Diagram 
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3000 K 


2700 K 


0.26 0.28 0.30 0.32 


( e r 7-step 

MacAdam ellipses 
(CFL) 


0.34 0.36 0.38 0.40 0.42 
X 


Table 1 - Nominal CCT Categories 


Nominal CCT 1J 

Target CCT and 
tolerance (K) 

Target Duv and tolerance 

2700 K 

2725 

± 

145 

0.000 

± 

0.006 

3000 K 

3045 

± 

175 

0.000 

i 

0.006 

3500 K 

3465 

± 

245 

0.000 

± 

0.006 

4000 K 

3985 

± 

275 

0.001 

± 

0.006 

4500 K 

4503 

± 

243 

0.001 

± 

0.006 

5000 K 

5028 

± 

283 

0.002 

± 

0.006 

5700 K 

5665 

± 

355 

0.002 

± 

0.006 

6500 K 

6530 

± 

510 

0.003 

± 

0.006 

Flexible CCT 
[2700 - 6500 K) 

7 2 > 

± 

AT 3 ' 1 

D 

± 

0.006 
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CCT- Duv chart 


Duv 



2000 3000 4000 5000 6000 7000 8000 


(CCT in log scale) 
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Correlated Color Temperature (CCT) 


Temperature [K] of a Planckian radiator whose chromaticity is 
closest to that of a given stimulus on the CIE ( u',2/3 v') coordinate. 


CIE 1960 u-v Chromaticity Diagram* 



u 


(CIE 15:2004) 

CCT is based on the CIE 
1960 ( u , v) diagram, which 
is now obsolete. 

CCT is valid within distance 
0.05 from the Planckian 
locus on the ( u ' , 2/3 - V ) 
diagram. (CIE 15: 2004) 
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CIE 15:2004 Colorimetry, 3rd Edition 


APPENDIX E, INFORMATION ON THE USE OF PLANCK’S EQUATION FOR STANDARD 
AIR 

According to the Planck’s law, the spectral radiance of a blackbody at thermodynamic 
temperature T[K] in a medium having index of refraction n is given by 

(E.1) 

where Ci = 2Khc 2 T c 2 = hc/k, h is Planck’s constant, c is the speed of light in vacuum, k is the 


t-ej (X, T} = — 


_ „ -2 i-5 

C-iD a 


71 


exp 


r 




nXT 


-1 


T should follow the current International Temperature Scale (ITS-90), therefore, 
c 2 — 1,4388 x 10' 2 m K. 


Therefore, in the current recommendation in CIE 15:200 4. colour tempera 
correlated colour temperature are calculated using Equ. E.1 with 
change from the previous practice. This recommendation may be sub 


n = 1 (exactly 1 ), 


ure and 
thus no 


ect to change in the future. 




CORM 2011 


III^B 


NIST 9 


Robertson (1968) 


Computation of Correlated Color Temperature and Distribution 
Temperature, Journal of the Optical Society of America, 58-11, 1968 



Fig. 1 + Method of interpolation to find correlated 
color tempera to re* 


Tvcitl HI. Maximum errors oi computed values of cornsbied 
color icnpmLuie, Litcd on use of the 31 i$otecnpmture lincj 
listed in Tahk II. 
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Direct approach (1) to calculate CCT and Duv 


Triangular solution 



RB Locus 
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0.7&?2i 
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(1 ) Create a table of CCT vs 
distance d, to BB locus on (u,v) 
coodinate. 

(2) Find the closest point in the 
table. 

(3) Solve the triangle for the 
neighboring 2 points 

i2 


X = 


d l -d l + V 

U m-\ U m+] 1 


21 


D.WZ3Z 

D.M537 

0.00449 

Q.W5JS 

C i 
o.oflaas 
0.00994 
0.01 101 


(u ’ V >T m+| Tx 


( u > v )t 


m-1 


r = t . ,).y 

A,v = [+Sign ] (<£_, -X 2 ) 


1/2 



Use Planck' s equation and color matching functions at 1 nm interval. 
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CCT Error in Triangular Solution 


CCT error in Triangular soiiution (1% step table) 



Error increases when the point is far from Planckian locus. 
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Direct approach (2) to calculate CCT and Duv 


Parabolic solution 



773? 

7/3ES 

2/8V 

28)5 

7843 

2871 

7900 

7929 

2033 

2988 
■1018 
3048 
3078 
3109 
31 40 
3172 


0.70109 
CL73EJES7 
0.238/6 
0.25761 
0.75847 
0.2 3335 
0.75 475 
0.7bZ31 S 
0.2320/ 
0.23100 
0.74993 
0.24890 
0.74788 
0.74088 
0.24383 
0.24486 


0.331 IS 
0.3 3 OHO 
0.33044 
0.33007 
0.34970 
0.34032 
0.34893 
0.34E3SS 
0.34813 
0.34775 
0.34735 
0.34691 
0.3 48 S3 
0.3481 1 
0.34369 
0.34326 


0.00 8 SO 
0.00S29 
0.004 1 1 
0.00297 
0.00193 
0.001 22 
0.00144 
0.00737 
0.0033/ 
0.00446 
0.00550 
0.00666 
0.00770 
0.000 85 
0.00994 
0.01 101 



CCT IK] 


(1 ) Create a table of CCT vs 
distance of, to BB locus on (u,v) 
coodinate. 

(2) Find the closest point in the 
table. 

(3) Parabolic fit for the neighboring 
3 points. 

d{I T ) = aT 2 +bT+C 

-b 

d(T)’ = 2aT +b = 0 v T 

' ' A X A 

2a 

D ay =[±sign](aT 2 +bT^+c) 
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CCT Error in Parabolic Solution 


CCT error in Parabolic solution (1 % step table) 



u 



o 

(J 



0.19 


0.21 


23 


0.25 


0.27 


Much better, but the problem is on or very close to Planckian locus. 
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Combined Solution 


Use Parabolic solution 
but, take the CCT of 
Triangular solution for 

D uv <0.002 


Error in Combined solution (1 % step table) 
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Most Accurate Version (cascade expansion) 


Calculation of CCT (Most Accurate version, 1 000 K - 20000 K) Y. Ohno 


INPUT 

X 

0.3127 


1 

0.3290 


O'. 1978 
0.3122 


1 5 % step table 


1000 

0.44301 

0.35462 

0.2537 

1150 

0.41559 

0.35725 

0.2224 

1 323 

0.38439 

0.35927 

0.1924 

1521 

0.35512 

0.36039 

0.1645 

1749 

0,32327 

0.36036 

0.1390 

201 1 

0.30412 

0.35893 

0.1 161 

2313 

0.28231 

0.35620 

0.0957 

2660 

0.26431 

0.35203 

0.0775 

3059 

0.24853 

0.34679 

0,0614 

351 8 

0.23527 

0.34060 

0,0470 

4046 

0.22430 

0.33332 

0.0342 

4652 

0.21535 

0.32677 

0.0223 

5350 

0.20313 

0.31972 

0.0127 

6153 

0.20237 

0.31292 

0.0046 

7076 

0.19731 

0.30655 

0,0057 

S 1 37 

0.19423 

0.30073 

0.0120 

9358 

0.19143 

0.29552 

0.0179 

10761 

0,18923 

0.29093 

0.0230 

123 75 

0.18751 

0.28694 

0.0273 

14232 

0.18615 

0.2B350 

0.0310 

16367 

0.18507 

0.28056 

0.0341 

1 8822 

0,18420 

0.27306 

0.0363 

21645 

0.18351 

0.27593 

0.0390 

24391 

0.18295 

0.27412 

0.0409 

23625 

0.18249 

0.27259 

0,0425 

32919 

0.1821 1 

0.27129 

0.0433 


match 

TC.m-1 ) 
T[m) 
T(m+1 ) 



0.1 5 % step table 


y 

V 

0.20032 

0.31 OB 6 

0.20077 

0.310B0 

0.20072 

0.31073 

0,20067 

0.31066 

0.20062 

0.31060 

0.20057 

0.3 106 3 

0.20053 

0.31047 

0.20048 0.3104D 

0.20043 0.31034 

0.20033 0.31027 

0,20034 

0.31021 

0.20029 

0.311014 

0.20024 

0.3100B 

0.20020 

0.31001 

0.20015 

0.309S5 

0.20010 

0.309BB 

0.20006 

0.309B2 

0.20001 

0. 309 7 6 

0.19997 

0.30969 

0.19992 

0.30963 


0.01 5 % step table 



u 

V 

Distance 

6494.3 

0.20043 

Q.31Q40 

0.003213 

6495,2 

0,20047 

0.31Q4Q 

0.003213 

6496.1 

0.20047 

0.31039 

0.003213 

6497.0 

0.20046 

0.3103B 

0.003213 

6498.0 

0.20046 

0.3103B 

0.003213 

6498.9 

0,20045 

0.31037 

0.003213 

6499.8 

0.20045 

0.31036 

0.003213 

6500.7 

0.20045 

0.31036 

0.003213 

6501.7 

0,20044 

0.31035 

0.003213 

6502.6 

0.20044 

0.31034 

0.003213 

6503.5 

0.20043 

0.31034 

0.003213 

6504.4 

0.20043 

0.31033 

0.003213 

6505.3 

0,20042 

0.31032 

0.003213 

6506.3 

0.20042 

0.31032 

0.003213 

6507.2 

0,20041 

0 . 3 ' 0 3 ' 

0.003213 

6508.1 

0.20041 

0.31030 

0.003213 

6509.0 

0.20040 

3.3' C3 = 

0.003213 

6510.0 

0.20040 

0.31029 

0.003213 

6510.9 

0.20039 

0.3102B 

0.003213 

651 1.8 

0.20039 

0.3 102 B 

0.003213 

6512.7 

0.20038 

0.31027 

0.003214 


match 


£512.7 Q.Q0321 

Parabolic solution 


TC.m-1} 

T(m) 

T(m+1) 


6503.5 0.00321 

Parabolic solution 


a 1.336E-05 




a 

0.001336 

b -3.77E-05 




b 

-0,003772 

c 1~,B B6E-05 




c 

0.001336 

A TI223E-OB 




A 

1.22E-03 

b -aoooi 59 




B 

-0,000 1 59 

C 0,5202373 




C 

0,520744 

CCT 6503.05 




CCT 

6503.03 

duv 0.00321 




duv 

0.00321 


Used as the reference 
for accuracy 
verification. 





T(m-1 ) u 

0.20032 




T( m-1 ) u 

0.20048 




TC.m-1 )u 

0.20044 




Ttm+1 > u 

0,19933 




T(m4l ) u 

0.20038 




T(m+1 ) u 

0,20043 




T(m-1 ) v 

0.31086 




T(m-1 ) v 

0,31040 




T(m-1 ) v 

0.31035 




T(m4l ) v 

0,30956 




T(m4l ) v 

0,31027 




T(m+1 ) v 

0.31034 




d 

0.00 




d 

0.00 




d 

0.00 




X 

0,00 




X 

0.00 




s 

0.00 




CCT 

6503,70 




CCT 

6503.03 




CCT 

6503,05 




duv 

0.00322 




duv 

0.00321 




duv 

0.00321 




CCT 

6504.39 




CCT 

6503.05 




Final CCT 

6503.03 



















0.2004 

0.3103 




0.2004 

0.3103 




0.2004 

0.3103 




Sign of Duv 

1 




Sign of Duv 

1 




Sign of Duv 

1 




Duv 

0.00322 




Duv 

0.00321 




Final Duv 

0.00321 



















CORM 2011 


NIST 16 


Conversion from (CCT, Duv) back to (x, y) 


Input: CCT 7"(K) 
Duv D uv 


1 ) Calculate (« 0 , v 0 ) of the Planckian 
radiator at T(K). 

2) Calculate («,, u x ) of the Planckian 

radiator at T+ 7"(K). r=0.01 (K) 

3) Calculate 

d?v = u t - u 0 
dv = V, - v 0 



u = u 0 + l) m - sin 6 

: u' = u 

= *'() + A, v- dv/idu 2 + dv 2 v' =l.5v 

v = v n + Av ' cos ^ a- = 9 1/'/(6 1/ - 1 6v' + 1 2) 

| 

= u 0 + /) uv - du/ydu 2 + dv 2 y = 2 v'/(3i/ - 8v' + 6) 


(Included in Revision draft of C78.377) 
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Accuracy of Most Accurate Version (4 stage) 
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Simple calculation from (x,y) or (u' ,v' ) to Duv 


Duv is normally calculated in the process of calculating CCT. 

Below is a simple approximation formula, without calculation of CCT. 


CIE 1960 u-v Diagram* 
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Simple calculation from (x,y) or (u' ,v' ) to Duv 


Accuracy of this method 


within 0.00001 in the range 
from 2600 K to 20000 K and 
Duv 0.000 ± 0.010 


within 0.0001 in the range 
from 2 1 60 K to 20000 K and 
Duv 0.000 ± 0.010 


> 

□ 

G 

e 

m 

L. 

e 

i_ 

LU 



0.00010 

0.00006 


0.00006 


0.00004 


0.00002 


- 0.010 

0.000 


0.010 


Duv 


0.00000 


CCT CK) 


(Included in Revision draft of C78.377) 
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Simple calculation from (x,y) or (u' ,v' ) to (CCT, Duv) 


06 19 CD u m 


A P = f( u - 0.292 ) 2 + (v -II.24) 2 

o ( =arclan((v -0.24 ) /(h — 0.292)), if c/, > 0, ;/ < 0, a = a x + n 

Am = kw, a + 

A iv =Ap “Ain 


Am ~~ ^06 + As <+ + A-i w 1 + A. 3 a + A2 + A i " + Ao 



For a < 2.54; 7j = 1/(A: 16 - a 6 + k l5 • a 5 + A: 14 - a 4 + & 13 - a 3 + & 12 - a 2 +k n - a + Ao) 

For a >2.54; T x = l/(Ar 26 • a 6 + As' + A4 ' fl4 + A 23 * fl3 + A2' q2 + Ai' a + Ao) 

For a < 2.54; A7^ = (A: 3 6' + As' a 5 + k 34 • a 4 + A3' + A2' ° 2 + Ai' a + Ao)* (Ab + 0.01)/ A * D uv /0.01 

For a >2.54; AT^j = 1/(A6' + As' + A 4 ' fl4 + A 3 ' « 3 + A 2 ' fl2 + Ai' a + Ao)* CAb +0.01)/ A * D uv /0.01 

T 2 = f - AT cl , c = log(r 2 ) 

For Duv >0; A7^ 2 = (Ao' ° 6 + As' + A4' c4 + A3' A + A2' c 2 + Ai' c + Ao) 

For Duv < 0; A T c2 = (k 66 - c 6 + k 65 • c 5 + k 64 • c 4 + k 63 - c 3 + k 62 - c 2 +k 6l - c + Ao)' |D UV / 0.03| 2 


Anal - A A A 


/ 

ki6 

ki5 

k/4 

ki3 

k/2 

k/1 

k/0 

0 

-3.7146000E-03 

5.6061400E-02 

-3.307009E-01 

9.750013E-01 

-1.5008606E+00 

1.115559E+00 

-1.77348E-01 

1 

-3.2325500E-05 

3.5700160E-04 

-1.589747E-03 

3.6196568E-03 

-4.3534788E-03 

2.1595434E-03 

5.308409E-04 

2 

-2.6653835E-03 

4.17781315E-02 

-2.73172022E-01 

9.53570888E-01 

-1.873907584E+00 

1. 96498025 1E+00 

-8.58308927E-01 

3 

-2.3524950E+01 

2.7183365E+02 

-1.1785121E+03 

2.51170136E+03 
6. 409763 56945E+0 

-2.7966888E+03 

1.49284136E+03 

-2.3275027E+02 

4 

- 1.731 364909 E+06 2. 748273293 5 E+ 07 

- 1. 8 1749963 507 E+08 

8 

-1.27141290956E+09 1.34488160614E+09 

-5.926850606E+08 

5 

-9.4353083E+02 

2.10468274E+04 

-1.9500061E+05 

9.60532935E+05 

-2.65299138E+06 

3.89561742E+06 

-2.3758158E+06 

6 

5.0857956E+02 

-1.321007E+04 

1.4101538E+05 

-7. 93406005 E+05 

2.48526954E+06 

-4.11436958E+06 

2.8151771E+06 
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Simple calculation from (x,y) or (u' ,v' ) to (CCT, Duv) 


Accuracy of this method 
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Conclusions 


• Practical calculation and conversion formulae for 
CCT and Duv have been developed. 

• Accuracies of some of the formulae will be further 
improved. 

• The use of CCT and Duv (rather than x, y or u’, v’ 
chromaticity coordinates) is recommended to 
specify the chromaticity of lighting sources. 
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THANK YOU for your 
attention. 


Contact: ohno@nist.gov 
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Proposed revision of ANSI C78.377 


Current version 



X 


Revision proposal 2 



X 


• All center points to be moved onto the Planckian locus. 

• This proposal is pending due to a need for vision 
experiments. 

• Anecdotes say people prefer below the Planckian locus. 

• NIST is funded by DOE to conduct vision experiments 
using STLF. 
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CCT (K) difference between 3 rd stage and 5 th stage 
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Summary 


• Duv is important for color quality of light sources. 

• Duv is often neglected in specifications. 

• Parabolic and triangle combined solution works well 
for CCT calculation. 

• 1 % step table provides enough accuracy 

(<1 K for 1000 to 10000 K, <2 K up to 20000 K, Duv± 0.03) 

• Most Accurate Version (cascade expansion), 

• Conversion from (CCT, Duv) back to (x,y), 

• Simple calculation from (x,y) or ( u ',v') to Duv, 

• Simple calculation from (x,y) or (u',v') to (CCT, Duv) 
have been developed. 






CORM 2011 


NIST 28 


